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Abstract Theory and observations favor stable helium burning as the most important means to produce fuel for superbursts on 
neutron star surfaces. However, all known superbursters exhibit unstable burning as well. This ambiguity prompted us to search 
for superbursts in data from the BeppoSAX Wide Field Cameras of ten luminous LMXBs, most of which do not exhibit normal 
type-I X-ray bursts. We found no superbursts and determine a lower limit on the recurrence time which varies between 30 and 
76 days (90% confidence). All recurrence time limits except one are longer than the observed recurrence time for GX 17+2. 
This difference can be understood if the mass accretion rate in GX 17+2 is several tens of percent higher than in the other 
sources; alternatively, the accreted material in GX-17+2 might be hydrogen deficient, leading to larger carbon yields than in 
the other sources. We compare our results to the latest models of superbursts. As our search method is indiscriminate of the 
burst ignition scenario, the recurrence time limits may also be applied to other bursts of similar duration and brightness. 
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1. Introduction 

Neutron stars in low-mass X-ray binaries (LMXBs) accrete 
hydrogen- and helium-rich matter from a sub-solar mass com- 
panion star. The gravitational energy that is released when the 
matter travels to the surface of the neutron star accounts for 
most of the observable X-ray radiation. The accreted matter 
forms a layer on the neutron star. When this layer is suffi- 
ciently compressed and heated, thermonuclear burning is ig- 
nited. The burning can be stable or unstable, depending on 
e.g. the mass accretion rate. Unstable nuclear burning re sults 
in a type I X-ray b urst (Grindlav et al. 1976; Belian et al. 1976; 
IWooslev & TaanJll976tlMaraschi & Cavalierelll977l) . A burst 
is observed as a fast (1-10 seconds) rise in the flux and a 
slow (10-1000 seconds), exponential-like decay. Many such 
bursts have been obse rved from about 8 sources (see e.g . 
lin 't Zand etai1l2004bl). For revie ws see iLewin et al] Jl993l) 
and lStrohmaver & Bildstenl 12003I) . 

Longer X-ray bursts have been observed that have a similar 
fast rise, but a much longer decay of several hours. These so- 
called superbursts release about a thousand times more energy 
than the normal type I X-ray bursts. So far, ten superbursters 
have been identified, with two exhibiting m ore than one su- 
perburst. For an observational overview, see iKuulkersI J2 004) 
(and see alsoliti 't Zand et all2004aIlRemillard et alJl2005l and 
lKuulkersl2005l for reports of new superbursters - ) . 

Superbursts are ascribed to unstable carbon burning (ini- 
tiated by 12 C + 12 C -> 20 Ne + a) in an 'ocean' of 



heavy elements which lies bet ween the freshly ac creted 
layer and the neutron star c rust JCumming & Bildstenl 120011 
IStrohmaver & Brownl [2002). The carbon and heavy ele- 
ments are thought to be the products of rp-process hy- 
drog en and helium burning in the accreted layer (see e.g. 
ISchatz et alJ 1 19991 120031) . This model has been fairly suc- 
cessful at reproduc ing the energetics and recurrence times of 
superbursts (Cumming & Bildsten 2001; Strohmaver & Browr 



own 



200i iBrowrJkoW lcooDer & Narava ^X^f lcumming et al 
2005), as well as the observed lightcurves and the quenching 
of norma l type I X-ray bursts for about a mo nth followinga su- 
perburst JCumming & Macbethl2004ICumming et alJ2005l) . 

An open issue is how to produce enough carbon to power 
a superburst. Carbon fractions of more than 10% are required 
to reproduce the observed lightcurves of superbursts and to 
achieve unstable ignition at the accretion rates of approxi- 
mately O.lMEdd inferred for most of the superbursters (Medd 
is the Eddington accretion rate). Calculations of rp-process 
burning show that stable rather than unstable H/He burning is 
required for such large carbon fractions JSchatz et alJ 12003). 
The reason is that carbon production is only possible by he- 
lium burning after the hydrogen has been consumed, since hy- 
drogen readily captures on carbon, processing it to heavy el- 
ements. Helium is rapidly consumed at the high temperatures 
reached during unstable burning. However, in stable burning, 
the helium burns at a reduced rate so that when the hydrogen 
has been completely consumed by rp-process burning, there is 
still some helium left that then burns into carbon. In this way, 
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there is an anti-correlation between the amount of carbon that 
is produced and the length of the rp-p rocess, or equivalen tly the 
average mass of the heavy elements (Sc hatz et al.l 2003). 

Observations support the prese nce of stable burning in su- 
perburst sources. In 't Zand et al. (2003) showed that super- 
burst sources preferentially have large values of alpha, the ratio 
of the persistent fluence between type I X-ray bursts to the flu- 
ence in the type I X-ray burst. Large values of alpha indicate 
that not all the nuclear fuel is consumed in burst s, and that sta- 
ble burning must occur between them. Corneli sse et alJ (|2003) 
studied BeppoSAX observations of nine type I X-ray bursters 
and found that the burst rate drops for luminosities exceeding 
approximately 10% of the Eddington limit, implying that sta- 
ble burning is occurring. This critical luminosity is close to 
the lowest luminosity of superburst sources, and may explain 
why superbursts have not been seen at accretion rates less than 
~ O.U/Edd. 

In contrast, the role of unstable burning and the 
heavy elements produced by the rp-process remains unclear. 
Observationall y, it is true that all supe rburst sources show type 
I X-ray bursts. ICumming & Bildsten (2001) pointed out that 
for a fixed heat flux emerging from the neutron star crust, the 
heavy nuclei reduce the thermal conductivity of the accumu- 
lating ocean, leading to higher temperatures and earlier carbon 
ignition tha n previously found in m odels of p ure carbon lay- 
ers (see e.g. Wooslev & Taam Il976l) . However. iBrownl d2004 
showed by considering the thermal state of the entire neutron 
star that in fact the thermal state of the core and crust of the 
neutron star is most important for setting the ignition column 1 . 
Because the flux emerging from the crust changes with the 
heavy element composition, the overall effect of the heavy ele- 
ments is much smaller than calculated bv lCumming & Bildstenl 
l l200ll) . Nonetheless, the heavy elements do have some effect 
on the ignition depth, and so it is possible that unstable burning 
could contribute to the presence of superbursts. 

One way to address this issue is to study the several persis- 
tently bright LMXBs which either do not or very rarely show 
type I X-ray bursts. The rapidly accreting source GX 17+2 
has already been important f or testing carbon ignition models. 
ICumming & Bildstenl (l200lh showed that superbursts should 
occur at rapid accretion rates near Eddington, but with re- 
duced energies and recurrence times. In 't Zand et al. (12004a) 
identified four superbursts from GX 17+2 which well-matched 
the predicted recurrence times, although their durations were 
longer than expected. In this paper, we report on a system- 
atic search for superbursts in BeppoSAX Wide Field Camera 
data of ten LMXBs with sufficiently high accretion levels. 
This includes the so-called "GX" sources: very bright LMXBs 
with luminosities close to the Eddington limit. We include six 
sources that never exhibited an X-ray burst as well as four X- 
ray bursters for comparison purposes, including the known su- 
perburster GX 17+2. We compare the behavior of the four most 
luminous sources Sco X-l, GX 340+0, GX 5-1 and Cyg X-2 
to the equally luminous superburster GX 17+2. For all these 



This is a very exciting result because it opens up the possibility 
of usin g superbu rst observations as a probe of the neutron star interior 
jBrownboO^ICooner & Naravarl2005l : ICumming et all2005t) . 



sources we find a constraining lower limit on the superburst 
recurrence time. 

2. Observations 

We search for s uperbursts in a rchival data from the Wide Field 
Camer as (WFC:I.Ta ger et alJl997l) on BeppoSAX dBoella et alJ 
Il997al) of ten X-ray sources. The two WFCs were identical 
coded mask telescopes pointing in opposite directions. They 
had 40° x 40° fields of view and a 5' angular resolution in the 
2-25 keV bandpass. They were active between 1996 and 2002. 
The WFCs carried out a program of semi-y early ca mpaigns 
on the G alactic center region. For a review see lin 't Zand et all 
(2004b). We consider nine sources near the Galactic center and 
Cyg X-2 (see Table [Q. The selected sources accrete mass at 
a comparable rate to the known superbursters. Some of the 
sources exhibit normal bursts while others do not. 

For each source there is 55 to 153 days of observation time 
accumulated over a period of six years (see Table [Q. Note 
that this is not the net exposure time, but includes data gaps 
due to Earth occultations (typically 36 minutes long and recur- 
ring every 96 minutes) and passages through the South Atlantic 
Anomaly (13 to 26 minutes long). During these data gaps there 
is no observation of the source. However, the duration of a data 
gap is short with respect to the expected duration of a super- 
burst. These data gaps are, therefore, not very detrimental to 
the detection of superbursts. Except for these short datagaps, 
an uninterrupted period of exposure time is typically 1.5 days 
long. 

To determine the persistent flux in the 0.1-200 keV range 
of GX 9+9 we use data from the BeppoSAX Narrow-Field 
Instruments (NFI). Of these instruments we use data from the 
Low Energy Concentrator Spectrometer (LECS: Iparmar etal] 
119971 0.1-10 keV), the Medium Energy Concentrator 
Spectrometers two and three (MECS; Boel la et alJ Il997bt 
1.3—10 keV) and the Phoswich Detection System (PDS; 
iFrontera et alJll997t 13-200 keV). For the other sources we 
obtain a measure of the flux from the literature (see Section[3}. 

3. Search for superbursts 

For each source we have 55 to 153 days of observation time 
while a superburst is typically observable for several hours. 
From the superbursts thus far detected we know that they are 
rare, so we expect to find at most a few new superbursts. The 
search for these consists of an automatic peak-search algorithm 
and a visual inspection of the data. 

We use a lightcurve in the full bandpass at three differ- 
ent time resolutions. The highest time resolution is 16 s. This 
resolution allows to discriminate between fast and slow rising 
peaks. A higher time resolution is not useful because the uncer- 
tainty in each lightcurve point becomes too large. The next time 
resolution is 0.02 days (about 0.5 hr). We find this resolution 
good to resolve superbursts at optimum sensitivity since short 
features such as normal bursts are averaged out while super- 
bursts are not. The lowest time resolution is 0.1 d. We employ 
it to measure the non-flaring persistent flux. At this resolution 
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Table 1. Superburst recurrence time of ten sources. For each source we give the time in days between the first and the last observation 
t S pan, the total observation time t b s including data gaps, the flaring fraction, the number of observed superbursts (SBs), whether a source 
exhibits normal bursts and the (lower limit of the) superburst recurrence time t lacul . 



Source tspan[d] t bs[d] flaring fraction # SBs burster t roC ur[d] 



Sco X-l 


1 OTA A"5 

18/0.03 


55.44 


16% 





no 


>30 










n 
\t 


no 


^>74 


GX 349+2 


2062.82 


129.71 


12% 





no 


>66 


GX9+9 


2062.57 


139.27 


8% 





no 


>76 


GX 354-0 


2062.79 


100.97 


10% 





yes 


>52 


GX5-1 


2062.79 


109.95 


11% 





no 


>55 


GX9+1 


2062.79 


103.45 


9% 





no 


>53 


GX 13+1 


2062.79 


104.30 


9% 





yes 


>55 


GX 17+2 


2062.80 


127.59 


16% 


4 


yes 


30 ± 15 a 


Cyg X-2 


1999.10 


109.67 


11% 





yes 


>60 



a Fromlin 't Zand et alj J2004ai. 



superbursts are largely averaged out while typical variations in 
the persistent flux are well sampled. 

The search algorithm proceeds as follows. At 0.02 d time 
resolution we check each data point whether it exceeds the per- 
sistent flux significantly. Each data point is assigned a mea- 
sure of the persistent flux through a linear interpolation of the 
preceding and following 0.1 d-resolution data point. We then 
check whether the 0.02 d-resolution data point exceeds the per- 
sistent flux at 4.4(7 significance. We choose 4.4a since at this 
level of significance we detect all superbursts of GX 17+2 
while limiting the number of false peaks. There are typically 
6000 points in a 0.02 d-resolution lightcurve. We expect the 
number of false peaks exceeding 4.4er due to Gaussian noise to 
be 0.06 . 

Using the lightcurve at 16 s time resolution we check these 
peaks visually for superburst characteristics: a fast rise and a 
slow, exponential-like decay. 

The lightcurve may contain featur es that have the charac - 
teristics of a superburst but are flares din 't Zand et alJl20 04a). 
Flares are thought to result from quick changes in the mass 
accretion rate. They are generally not expected to be iso- 
lated events, but to come in multiples during flaring episodes. 
Therefore we exclude obvious flaring periods. We identify flar- 
ing episodes visually from the light curve. The fraction of the 
lightcurve that shows flaring is presented in Tabled F° r the 
identific ation of a flaring episode we use a similar argument as 
lin 't Zand et alJ (|2004a): an interval between two consecutive 
data points at 0.1 d-resolution is identified as a flaring episode 
when the root mean square of the count rate at 0.02 d time res- 
olution in that interval exceeds the average root mean square 
of all intervals by a factor of two. The flaring fractions are on 
average 11%. 

Using this method we reproduce the four superbursts found 
previously for GX 17+ 2 as well as the two p eaks that were 
determined to be flares Jin 't Zand et alJ Eo04a). For the other 
sources we find no superbursts and we determine a lower limit 
on the superburst recurrence time t lecm . 



4. Recurrence time 

Thus far only two out of ten superbursters exhibited more than 
one superburst. These are therefore the only sources for which 
a recurrence time has been directly obtained. From GX 17+2 
four superbursts have been observed with an average recur- 
rence time of 30 ± 15 days and from 4U 1636-536 three su- 
perbu rst were observed wit h recurrence times of 1.75 and 2.9 
years (Kuulkers et a ll2004 . 

To determine the lower limit of the recurrence time for the 
other nine sources here investigated, we perform Monte Carlo 
simulations of series of superbursts and cross check those with 
the WFC observation schedule. We generate the series of on- 
set times, for a given t ICCUI , by first randomly picking the time 
of the initial superburst between and £ recU r- Subsequent su- 
perburst times are sampled from a Gaussian distribution with 
a standard deviation a, centered around a time t lccuI after the 
previous superburst. We use the Gaussian distribution to model 
the variability in the recurrence time resulting from the vari- 
ability of the mass accretion rate as measured through the flux. 
For the nine sources we find from the normalized root mean 
squared of the flux that a is between 0.06 and 0.3 i raC ur- 

The superburst exponential decay times are randomly sam- 
pled between 1 and 6 hours (see review by Kuulkers 2004). As 
we will show in more detail below, we can take the decay time 
independently from the wait time since the previous superburst, 
because we consider decay times longer than the duration of the 
data gaps. We assume the peak flux to be equivalent to half the 
Eddington limit of each source. This implies that any super- 
burst is detectable with the WFC for as long as the e-folding 
decay time. 

A given superburst is considered to be detected if at any 
time of its e-folding decay duration the instrument was observ- 
ing the source according to the WFC observation schedule, tak- 
ing into account all data gaps, including those due to Earth oc- 
cultations and passages through the South Atlantic Anomaly, 
as well as flaring episodes. 

For each source and i rocur we perform 100,000 Monte 
Carlo simulations. With this number of simulations one expects 
an uncertainty of 0.3% in the obtained value. We determine the 
percentage of simulations in which at least one superburst is 
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Length superburst (hours) 

Figure 1. Lower limit to the recurrence time at 90% confidence 
as a function of assumed superburst duration in Monte Carlo 
simulations. The lower limit drops steeply when the burst du- 
ration is of the order of the length of the data gaps or less. Note 
that our search is insensitive to bursts shorter than 0.5 hour due 
to our choice of lightcurve time resolution (see Section^. 



detected. We define the lower limit to the recurrence time to 
be that value of t recur for which there is a 90% probability of 
detecting a superburst. The results are provided in Tabled We 
find lower limits for t tecuI between 30 and 76 days. 

In Figure H we illustrate the dependence of the superburst 
detection on the duration from the Monte Carlo simulations. 
For this we extend the duration range to below 1 hour. This 
clearly shows that for superburst durations above 1 hour, the 
presence of short datagaps is not detrimental to the detection of 
superbursts. 

Our aim is to look for superbursts, but in fact our search 
is not limited to carbon flashes. The search method we use 
is indiscriminate of the ignition scenario. Other X-ray bursts 
of sufficient duration and brightness will also be found using 
our method. Long duration bursts have b een observed from 
severa l sources (see [gumming e t alJl2005l and iKuulkers et alJ 
120021) . However, the exponential decay time of these bursts is 
< 600 s, which is shorter than the lightcurve time resolution 
of 0.02d that we use (see Section 0. Therefore we would not 
have been able to detect these bursts. Note that if a lightcurve 
at a shorter time resolution is used, such that these long burst 
are detectable, the limits that can be placed on the recurrence 
time are far less constraining than for superbursts, since the 
burst duration is in this case shorter than the data gaps in our 
observations (see Figure [Q. Nevertheless, the constraints we 
find on the recurrence time apply not only to superbursts, but 
to any type of thermonuclear flashes which are long and bright 
enough to be detected. 

5. Mass accretion rate 

The unstably burning carbon layer to which superbursts are as- 
cribed is formed by the burning of accreted hydrogen and he- 
lium. The rate at which the hydrogen and helium are accreted 



from the companion star plays an important role in whether 
they burn stably or unstably. 

A measurement of the mass accretion rate is most directly 
done through the X-ray luminosity because the power for the 
radiation is provided by the liberation of gravitational energy of 
the accreted matter. However, this is not unambiguous because 
it is somewhat uncertain what fraction of the liberated energy 
goes into radiation. It is generally thought that some fraction 
goes into the kinetic energy of jets. Nevertheless, in the case 
of neutron star accretors it is not unreasonable to assume that 
most of the energy goes into radiation (see e.g. lFendeil2002h . 

Another complication involves the determination of the lu- 
minosity: distances to LMXBs — needed to translate fluxes to 
luminosities — are generally poorly known (except for systems 
in globular clusters). This applies particularly to many of the 
systems that we study in this paper: persistent LMXBs outside 
globular clusters without type-I X-ray burst behavior. Systems 
that exhibit bursts with Eddington-limited fluxes, as diagnosed 
through photospheric radius expansi on, provide a means to find 
distances with « 30% accuracy ("cf.. lKuulkers et alT20 03i. and 
transients that turn quiescent may reveal the optical light of the 
companion star without contamination by optical light from the 
accretion disk. Neither of these two conditions apply to most of 
our systems. 

A third complication is that for many systems we do not 
have a broad-band measurement of the X-ray spectrum that 
is averaged on timescales longer than that of the typical vari- 
ability to enable a good bolometric correction. Broadband 
BeppoSAX X-ray spectra have been taken in many cases but 
almost never for long enough. This situation may change soon 
with the results of the deep exposures of INTEGRAL on the 
Galactic bulge. 

We attempted to obtain the best broadband flux measure- 
ments for the ten systems here studied, see Table [2] For nine 
sources this pertains to values found in the literature from 
RXTE, BeppoSAX and INTEGRAL observations, and for one 
system, GX 9+9, we perform the analysis ourselves from 
BeppoSAX data. The data were taken on April 8, 2000, with 
exposure times of 23, 50 and 22 ks for LECS, MECS and PDS, 
respectively. Background spectra are taken from blank field ob- 
servations for LECS and MECS, at the same detector positions, 
and from off-source pointings with the PDS. We rebin the spec- 
tra to limit the oversampling of the spectral resolution to 40% 
of the full- width at half maximum, and to obtain at least 15 
photons in each spectral bin. During the spectral analysis, a 
systematic error of 1% is added in quadrature to the statisti- 
cal error per spectral bin and the normalization between the 
three instruments is allowed to vary. We find the generic LMXB 
spectral model to provide g ood fits to the data. This model 
(see e.g. ISidoli etalJ EoOl) consists of a multi-temperature 
disk black body jMitsuda etai]ll984t iMakishima et aljfl 986) 
in combination with a c o mptonized spectrum (jj itarchuk 1994; 
lHua&Titarchuklll995t iTitarchuk & Lvubarskiil 1 1995b . both 
absorbed by cold i nterstellar matter following the model by 
M orrison & McCammon lll983l) . The results of the fit are pro- 
vided in Table |3] 

The flux measurements are based on observations which 
were typically performed during 1.3 days. From the one-day 
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Table 2. Flux and luminosity estimates of the 10 investigated sources. We give the average persistent flux F in the energy range 
indicated below, the relative variability of the flux as determined from ASM/RXTE measurements, distance d, luminosity L in units of the 
Eddington luminosity LEdd — 2 • 10 38 erg s _1 and the ignition column depth y. We indicate the uncertainty when given, otherwise we assume 



a 2% uncertainty. The distances are taken from| 


Christian & Swank 


1 19971) unless indicated otherwise. Christian & Swank take the uncertaintv 


in the distance to be 30%.I.Tonker & Nelemansl 


(2004 


•) indicate that the uncertainty in the distances they determined is "probably very large". 


Therefore we use the 30% from Christian & Swank. The uncertainty in the luminosity is dominated by the uncertainty in the distance. Since 



Lad 2 , the uncertainty in L is 60%. As the ignition column depth is proportional to the luminosity (see equation^, it has the same uncertainty 
of 60%. In the case of GX 17+2 this has to be added to the indicated uncertainty in y from the spread in the observed recurrence time. 



Source 


F [lO^ergs^cm" 2 ] 


variability 


d [kpc] 


L/L-Edd 


y [lO^gcm" 2 ] 


Sco X-l 




13% 


2.8±0.3 a 


0.9±0.3 b 


>1.8 


GX 340+0 


2.67 c 


15% 


<11 


<1.9 


>9.4 


GX 349+2 


2.45 d 


16% 


5.0 


0.4 


>1.8 


GX 9+9 


0.8 ± 0.012 6 


14% 


5.0 


0.1 


>0.6 


GX 354-0 


0.037 f 


37% 


4.5 g 


0.05 


>0.2 


GX 5-1 


2.56±0.4 h 


8% 


<9.0 


<1.2 


>4.4 


GX9+1 


2.F 


13% 


5 1 


0.3 


>1.0 


GX 13+1 


0.803 ± 0.004 j 


7% 


7.0 


0.2 


>0.7 


GX 17+2 


2.26 k 


12% 


7.5 


0.8 


1.6 ±0.8 


Cyg X-2 


1.3 1 


23% 


8.0 


0.5 


>2.0 



a Dista nce from parallax meas urements I Bradsha w et alll999li 
b From Bradsh aw et all <2003l) (2.0-18.2 keV) 
c From lLavagetto et al J |2004t) (0. 1-200 keV) 
d From llaria et alJfcOoA CO. 1-200 keV) 

c Fro m observations with Bep poSAX NFI (0.1-200 keV). See Tabled 
f From lGalloway et alJl2003l) (2-60 keV) 

g Dista nce fromljon ker & N elemansl 120041) for hydrogen accreting neutron star. 

h From Paizis et al. (2003) (1.5-40 keV). Note that the flux in the 12-20 keV range is obtained through linear interpolation. 
' From Iaria et alj <2005t) (0.12-18 keV) 

J From lCorbetH2003li CI 5-12 keV) and lRevnivtsev et alJ <2004l) (18-60 keV). Flux in 12-18 keV range obtained through linear interpolation. 

Total e nergy range: 1.5-60 keV . 

k From in 't Zand et alJ <2004ah (0.1-200 keV) 

1 From lDi Salvo et al] fc002) (0.1-200 keV) 



Table 3. Flux of GX 9+9 from BeppoSAX NFI observations. 
As model we use an absorbed disk black body in combination with a 
comptonized spectrum. We provide the parameters from the best fit of 
the model to the observed spectrum and the unabsorbed 0. 1-200 keV 
flux. 





(0.231 ± 0.003) 10 22 crn^ 2 


kTsB 


0.853 ± 0.002 keV 


FIuxbb (0. 1-200 keV) 


(0.33 ± 0.014) 10" 8 ergs" 1 cm" 2 


kT 


1.055 ± 0.003 keV 


kT e 


2.458 ± 0.004 keV 


T 


14.10 ±0.06 


Fluxcomptt (0. 1-200 keV) 


(0.50 ± 0.04) 10~ 8 ergs" 1 cm" 2 


X 2 /dcf 


166/147 


Flux (0. 1-200 keV) 


(0.83 ± 0.013) 10" 8 ergs" 1 cm -2 



averaged ASM/RXTE lightcurves we determine the fractional 
rms variability for each source, see Table |2] We applied the 
best distance estimates from the literature and determined the 
ratio of the resulting luminosity and the Eddington luminos- 
ity of a canonical neutron star (i.e. iEdd = 2 • 10 38 ergs _1 



for a neutron star with a mass of 1.4M , a 10 km radius 
and a hydrogen-rich photosphere). Due to the three sources 
of uncertainty eluded to above, the resulting ratios are never 
very accurate, and we are only able to make a crude distinc- 
tion between the sources close to one tenth and 100% of the 
Eddington limit. This inference is in line with the fact that 
the five sources close to the Eddington limit are the only ones 
that trace out a Z-shaped curve in a har dness-intensity dia- 
gram (e.g. JHasinger & van der Klislll9 89). They are Sco X-l, 
GX 340+0, GX 5-1, GX 17+2 and Cyg X-2. The others are 
so-called Atoll sources. 



6. Discussion 

We analyzed a large volume of data from the Wide Field 
Cameras on BeppoSAX of ten luminous LMXBs to search for 
superburs ts. Except for tho se already reported from GX 17+2 
bv lin 't Zand et al.l ll2004al) . none were found and we set lower 
limits on the recurrence time of superbursts that range from 30 
to 76 days. 

One of the conditions for a superburst to occur is a suffi- 
ciently large column depth of accumulated matter. The ignition 
column depth y for a 1.4 M© neutron star with a 10 km ra- 
dius, assuming a gravitational redshift z = 0.31, is given by 
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Figure 2. Superburst recurrence time as a function of the mass ac- 
cretion rate in units of the Eddington limited mass accretion rate. We 
show models from Cum ming et alJ 120051) and observed recurrence 
times including our lower limits. The solid curves show results for 
a model with a disordered crust, a composition of 56 Fe and a 20% 
carbon abundance, without Cooper pairing for four different core neu- 
trino emissivities. The dashed curve is for a model with heavier com- 
position (A — 104), the dot-dashed curve has a higher crust conduc- 
tivity and the dotted curve includes Cooper pair neutrino emission. 
The curves end at lower m where stable carbon burning begins. 
Note that the uncertainty in the mass accretion rate is at least 60% for 
most sources (see Table|2j- 

Jin 't Zand et al.l2004al) 

y = 2.0.10 11 gcm- 2 f-|^) (1) 
V 30 days/ \m E dd/ 

with m(Edd) the (Eddington limited) mass accretion rate. At a 
given accretion rate, the lower limit to the recurrence time gives 
a lower limit to the ignition column depth. Using Equation Q 
we calculate y (see Table |2). 

ICumming et alJ l2005t I investigate several carbon ignition 
models for superbursts. From these the recurrence time as a 
function of the mass accretion rate is determined. In Figure |2] 
we compare our lower limits to t ICCUI as well as the observed 
recurrence times of 4U 1636-536 and GX 17+2 to these mod- 
els. Looking at what theory predicts, the constraints that we 
can put on the superburst recurrence time are only meaning- 
ful for sources that accrete close to the Eddington limit. For the 
sources that accrete close to 10% of the Eddington limit the the- 
ory predicts for all models a higher t lecnl than our lower limit. 
Sco X-l, GX 340+0, GX 5-1 and Cyg X-2 accrete matter at a 
rate comparable to GX 17+2, namely close to the Eddington 
limit 

Cumming et all d2005h discuss the dependence of the ig- 
nition of carbon on the heat flux into the superbursting layer 
from the deeper lying crust: Q/,. At a given accretion rate, 
the recurrence time is lower for higher Qi,. Therefore a lower 
limit to irccm, or alternatively to y, gi ves an upper limit to Q b 
for a given model (see Figure 20 in ICumming et alJ (Eo05)). 
Comparing the lower limit to the ignition column depth of the 



five brightest sources to the results of ICumming et alJ d2005l) 
we find Qb < 0.1 MeV per nucleon for iron composition. 
However, the large uncertainty in the accretion rate due to 
the large uncertainty in the distance to these sources pre- 
vents us from placing stringent constraints on the models of 
ICumming et al.ld2005l) . 

We found the lower limit to the superburst recurrence time 
for GX 340+0, GX 5-1 and Cyg X-2 to be about two times 
the recurrence time of GX 17+2. One of our goals was to learn 
more about the dependence of superbursts on normal bursting 
behavior. However, we cannot explain the difference in t rccul 
by the difference in normal type I bursting behavior as both 
GX 17+2 and Cyg X-2 exhibit bursts, while GX 340+0 and 
GX5-1 do not. 

From Figure[2]we see that, depending on the model, a 35% 
difference in the mass accretion rate can lead to a two times 
longer recurrence time. The uncertainty in the mass accretion 
rate is 60%. Therefore the difference in (lower limit on) t Tecnt 
could be due to GX 17+2 accreting faster than GX 340+0, 
GX 5-1 and Cyg X-2. 

Another possible explanation for the difference in t rccul is 
that GX 17+2 may differ from the other sources in some aspect 
which would lower the theoretical curves for this source. The 
core temperature of GX 17+2 might be higher, which would 
lead to a lower ignition column depth and the refore a lowe r 
recurrence time (see e.g. Figure 11 in ICumming et alJ Eo05). 
Alternatively there may be a difference in the hydrogen and 
helium abundance of the accreted matter. A larger initial he- 
lium abundance will lead dir ectly to a l a rger carbon yield from 
rp-process burning (Schatz et al. 2003). Cumming & Bildsten 
(2001) find that the amount of carbon required to achieve un- 
stable ignition at accretion rates close to the Eddington limit 
is approximately 3-5%. The presence of observable super- 
bursts in GX 17+2 might indicate that the accreted material 
in this source is hydrogen deficient, allowing sufficient car- 
bon to be made, whereas too little carbon is made in other 
sources to result in a superburst. Another effect, which goes 
in the opposite direction, is that a lower initial hydrogen abun- 
dance will lead to a less extensive rp-process, giving less mas- 
sive heavy elements, which increases the carbon ignition depth 
because of the lower opacit y of the accumulati ng fue l layer 
dCumming & Bildstenll20oil) . However, iBrownl d2004h shows 
that this effect is smaller than initially thought. At present only 
one superburster, the ultracompact system 4U 1820-30, is as- 
sumed to accrete hydrogen-poor matter. In this c ase, the ex- 
treme ly short orbital period of only 1 1 .4 minutes dStella et alJ 
1987) implies that the hydrogen mass fraction is X < 10%, 
consistent with the very luminous and en ergetic superburst ob- 
served from this source dCummingl2003l) . Another case of a su- 
perburster accreting hydrogen-deficient matte r is possibly the 
ultrac ompact system 4U 0614+091 (see e.g. iNelemans et all 
2004) of which most recently a superburst was observed 
dKuulkersfeOoll) . 

7. Conclusion 

We searched for superbursts from ten bright LMXBs. None 
were found but the four known superbursts from GX 17+2. For 
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the other sources we obtain lower limits to the recurrence time 
and the ignition column density. Comparing this lower limit 
for the five most luminous sources, Sco X-l, G X 340+0, GX 5- 
1, GX 17+2 and Cyg X-2, to the models by Cum ming et alJ 
(2005]) Q b < 0.1 MeV per nucleon is required. Due to the large 
uncertainty in the distance to these sources, the uncertainty in 
the mass accretion rate is large as well. This pr events us f rom 
placin g firm constraints on the predictions bv lCumming et all 
(2005). As our search method is indiscriminate of the burst ig- 
nition scenario, the recurrence time limits may also be applied 
to other bursts of similar duration and brightness. 

The four luminous non-superbursters have a mass accre- 
tion rate comparable to the luminous superburster GX 17+2. 
For three of these sources we find a lower limit to the recur- 
rence time that is twice as long as the observed average < rccur 
of GX 17+2. Most likely this difference is due to GX 17+2 ac- 
creting faster than the other sources and therefore GX 17+2 has 
a lower t Iecm , However, we are unable to determine whether 
this is the case because of the large uncertainty in the mass ac- 
cretion rate. 

We found no superbursts in the six non-bursting sources we 
considered, which confirms the current observational view of 
superbursters, i.e. all known superbursters exhibit normal type 
I X-ray bursts as well. This supports the suggestion that both 
stable and unstable burning of hydrogen and helium are neces- 
sary to produce the fuel for superbursts. 

To better constrain models of neutron star interior physics, 
more long term monitoring observations of (candidate) super- 
bursters are necessary. As superbursts have a long recurrence 
time, a long exposure time is req uired for each obse rved source. 
The planned mission MIRAX jBraga et alJl2004 will be the 
most suited for such a task. MIRAX will have on board the 
spare flight unit of the BeppoSAX WFCs as well as two hard 
(10-200 keV) X-ray cameras. The aim of MIRAX is to per- 
form continuous broadband observations of the Galactic center 
during nine months per year. This will provide a high probabil- 
ity to observe any superbursts that have a recurrence time less 
than approximately these nine months, which includes the pre- 
dicted recurrence time of many models shown in Figure |2] for 
sources which accrete close to the Eddington limit. The latest 
observations of superbursts ( Remillard et al. 2005; Kuulkers 
l2005h were done with RXTE/ASM, which observes 80% of 
the sky each 90 minutes. Due to a time resolution of 90 s the 
short (1-10 s) burst rise cannot be observed, which makes it 
more difficult to distinguish between a superburst and a flare. 
Observations from MIRAX will have a time resolution better 
than 0.5 ms which is more than adequate to observe the rise of 
superbursts. 
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